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1 .  Introduction  and  characterizations 

A  stochastic  process  X  =  (X(t))tg^  is  self-similar  with  parameter  H  >  0 
(H-ss)  if  X(c*)  =  cH  X(*)  for  all  c  >  0,  and  has  stationary  increments  (si) 

if  X(*  +  b)  -  X(b)  =  X(*)  for  all  b  €  R,  where  =  means  the  equality  of  all 

finite-dimensional  distributions.  Every  H-ss  si  process  is  stochastically 
continuous  [VI],  A  real-valued  stochastic  process  X  is  a-stabie  if  all 
finite-dimensional  distributions  are  a-stable  where  0  <  a  <  2,  and  when 
a  =  2.  it  is  Gaussian. 

The  only  H-ss  si  Gaussian  process  is  the  fractional  Brownian  motion  (see 
(MV]).  On  the  other  hand,  when  0  <  a  <  2,  many  H-ss  si  a-stable  processes 

have  been  recognized  in  the  past  ten  years  (see  [KM]).  Among  these,  the 

linear  fractional  Levy  motion  seems  to  be  the  most  important  example  of  an 
H-ss  si  a-stable  process  for  0  <  a  <  2.  and  its  properties  have  been  well 
studied  (see,  e.g. ,  [CM],  [ST],  [ALT]). 

As  discussed  in  our  previous  paper  [CM],  to  any  stochastically 
continuous  si  a-stable  process  X  =  (X(t))t€j^  with  1  <  a  <  2  correspond  two 
stationary  a-stable  processes  =  (Y^(t))t€JR  =  (^R(t))t€R  in  the 

following  way: 

f  0  r  t 

(1.1)  YL(t)  =  eU[X(t)-X(t+u)]  du  =  X(t)  -  eU-t  X(u)  du.  t  €  R, 

'  —CO  j  —CO 


(1.2)  Y  (t)  =  j  e  U[X( t)-X( t+u)]  du  =  X(t)  -  f  e£  U  X(u)  du. 
J0  Jt 


t  €  R. 


A  straightforward  calculation  gives  for  all  s  <  t. 
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(1.3) 

X(t)  -  X(s)  =  YL(t)  -  YL(s)  + 

[  Y. (u)  du 

Js 

and 

(1.4) 

X(t)  -  X(s)  =  YR(t)  -  YR(s)  - 

f  YR(u)  du. 

((1.3)  is  given  in  (2.2)  on  p.  307  in  [CM],  where  the  minus  sign  before  the 
integral  should  be  replaced  by  a  plus  sign  as  in  (1.3)).  (1.4)  is  derived 

from  (1.2)  in  the  same  way  as  (1.3)  is  derived  from  (1.1)  in  [CM]. 

The  existence  of  the  integrals  in  (1.1)  and  (1.2),  and  the  derivation  of 
(1.3)  and  (1.4)  from  (1.1)  and  (1.2),  respectively,  are  justified  only  when 
1  <  a  <  2.  (but  both  are  expected  to  be  valid  for  0  <  a  <  1  as  well). 
Therefore,  whenever  or  YR  is  considered,  we  assume  1  <  a  <  2.  The 
correspondence  between  X  and  Y^  is  one-to-one  and  so  is  the  correspondence 
between  X  and  YR.  Furthermore,  it  follows  from  (1.1)  and  (1.3)  that  either 
one  of  X  and  Y^  is  expressed  linearly  in  terms  of  past  values  of  the  other, 
so  they  are  left-equivalent  in  the  sense  that  if  if^(X.t)  is  the  closure  in 
probability  of  all  finite  linear  combinations  of  the  values  to  the  left  of  t: 
(X(u),  u  £  t),  then 

^(X.t)  =  Sf^(Y^.t),  for  any  t  €  K. 

Similarly  it  follows  from  (1.2)  and  (1.4)  that  if  S£R(X,t)  is  the  closure  in 
probability  of  all  finite  linear  combinations  of  the  values  to  the  right  of 
t:  (X(u) ,  u  £  t),  then 

^(X.t)  =  ^(Yjj.t),  for  any  t  €  IR. 

We  will  call  Y^  and  YR  the  left-equivalent  and  right-equivalent  stationary 
process  of  the  si  process  X. 

One  of  the  most  important  classes  of  stationary  stable  processes 


consists  of  moving  averages: 
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(1.5) 


Y(t) 


I 


g(t-u) 


dM  (u). 

a  ' 


t  €  IR. 


where  0  <  a  i  2.  g  €  La(IR).  is  a-stable  Levy  motion  (Brownian  motion  when 
a  =  2).  i.e.,  has  stationary  independent  a-stable  increments  with 


exp{- | t-s | 0^/2} 


if  a  =  2. 


(1.6)  E[exp{i0[Ma(t)-Ma(s)]}]  = 


exp{- | t-s | |0 |a[l  -  i/3  sgn(0)  tan(ra/2)]j 

if  0<a<2,  a  /  1 , 


exp{-  { t-s  j  ( 0  | [  1  +  i/3(2/w)sgn(0)£n|0|]} 

if  a  =  1, 


where  |p|  £  1.  We  note  that  for  each  a  €  (0  2],  Mq  has  1/a-ss  increments. 

When  g(x)  =  0  for  x  <  0.  Y  is  a  nonanticipating  moving  average  process 

as  it  is  expressed  in  terms  of  the  past  increments  of  the  stable  L6vy  motion 

Mq:  Y(t)  =  g(t-u)  dMQ(u);  and  when  g(x)  =  0  for  x  >  0.  Y  is  fully 

anticipating  moving  average  as  it  is  expressed  in  terms  of  the  future 

00 

increments  of  the  stable  Levy  motion  M^:  Y(t)  =  Jt  g(t-u)  tffl^u). 

The  linear  fractional  Levy  motion  (Aq  ^(a.b;  t))tQp  is  an  H-ss  si 
a-stable  process  defined  for  0  <  H  <  1,  0  <  a  £  2,  H  /  1/a,  a.b  €  IR.  by 

(1.7)  Aa  H(a,b;t)  =  (a[( t-u)”“1/a  -  (-u)”"1/a] 

+  b[(t-u)”"1/a  -  (-u)”'1/a]}  dMa(u), 

with  the  convention  0  =0  even  for  -r  <  0,  where  we  assume  /3  =  0  in  (1.6)  if 

a  =  1.  When  a  =  2,  the  processes  Ag  ^(a.b;*)  are  multiples  of  the  fractional 
Brownian  motion  for  all  a  and  b.  However,  when  0  <  a  <  2,  different  lines 
through  the  origin  of  the  (a.b)  plane  give  different  processes  ([CM],  [ST]). 

We  now  assume  1  <  a  <  2.  As  shown  in  [CM],  the  left-equivalent 
stationary  a-stable  process  Y^(a,b;*)  of  the  linear  fractional  L6vy  motion 
A a  jj(a.b;*)  is  a  moving  average 
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f  °° 

(1.8) 

YL(a,b;t)  =  g&  fe(t-u)  dMQ(u).  t  €  IR. 

where 

,  »  H-l/a  -x  (*+  v  H-l/a  ,  ] 

ga ,b(x)  *  a  x+  “  e  J0  e  v  dv 

(1.9) 

H-l/a  -x  f°  -v  H-l/a  , 

+  b  x_  -  e  e  v  dv 

Jx_ 

=  ;  a  gl(x)  +  b  g^x). 

When  Y,  (a.b;*)  is  a  nonanticipating  moving  average  of  M  .  namely  g  ,  (x)  =  0 

for  x  <  0,  then  b  =  0,  since  for  x  <  0,  g^(x)  =  0  but  ggfx)  /  0.  The  linear 

fractional  Levy  motion  with  b  =  0,  A  „(a,0;*),  is  the  one  introduced  in  [TW] 

a,n  " 

and  is  called  left  linear  fractional  L6vy  motion.  Hence,  given  a  linear 
fractional  L6vy  motion,  if  its  left-equivalent  stationary  process  is  a 
nonanticipating  moving  average,  then  it  is  a  left  linear  fractional  Levy 
motion. 

In  this  paper  we  generalize  this  characterization,  replacing  the  class 
of  linear  fractional  Levy  motions,  by  the  class  of  all  H-ss  si  a-stable 
processes  whose  left-equivalent  stationary  stable  processes  are 
nonanticipating  moving  averages. 

A  similar  situation  prevails  for  the  right-equivalent  stationary  process 
YR.  As  in  (1.8)  and  (1.9),  the  right  equivalent  stationary  a-stable  process 
Yd  of  the  linear  fractional  Levy  motion  A  u  is  a  moving  average 

(1.10)  YR(a,b;t)  =  [  ha  b(t-u)  dMJu), 

*'  —00  ‘ 


where 
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(1.11) 


ha>b(x)  =  a 


H-l/a  x 


-  e 


L 


-v  H-l/a  , 
e  v  dv 


+  b 


f 

J0 


H-l/a  x  f  -  v  H-l/a  , 
x  -  e  |  e  v  dv 


=:  a  hj(x)  +  b  hgfx) 


When  YD(a,b;*)  is  a  fully  anticipating  moving  average  of  M  ,  namely 

K  Ct 

h  ,  (x)  =  0  for  x  >  0.  then  a  =  0,  since  for  x  >  0.  h0(x)  =  0  but  h,(x)  /  0. 
a,  b  2  1 

Thus,  given  a  linear  fractional  L6vy  motion,  if  its  right-equivalent 
stationary  process  is  fully  anticipating  moving  average,  then  it  is  a 
right  linear  fractional  L6vy  motion;  and  we  also  generalize  this 
characterization  to  general  H-ss  si  a-stable  process  as  follows. 

Theorem  1 .  Fix  l<a<2,  0  <  H  <  1,  H  /  1/a.  Let  X  be  a  nondegenerate  H-ss 
si  a-stable  process. 

(i)  If  its  left-equivalent  stationary  process  Y^  is  a  nonanticipating 
moving  average,  then 

X(*)  =  A  „(a,0;*)  for  some  a  /  0. 
a ,  n 

(ii)  If  its  right-equivalent  stationary  process  Y^  is  a  fully 
anticipating  moving  average,  then 

X(*)  =  ^(0,b;*)  for  some  b  /  0. 

We  suspect  that  the  linear  fractional  L6vy  motions  ^  of  (1.7)  are  the 
only  H-ss  si  a-stable  processes  whose  left-  (or  right-)  equivalent  stationay 
stable  processes  are  moving  averages.  However,  we  are  not  able  to  prove  this 
characterization  at  present. 

It  is  seen  from  (1.9)-(1.12)  that  the  left-equivalent  stationary  process 
Y^  of  a  nondegenerate  linear  fractional  L6vy  motion  cannot  be  a  fully 
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anticipating  moving  average,  because  if  g  ,  (x)  =  0  for  x  >  0,  then  a  =  b  =  0 

3 1  O 

and  g&  k  =  0;  and,  likewise,  the  right-equivalent  stationary  process  Y^ 
cannot  be  a  nonanticipating  moving  average. 

A  related  characterization  of  the  linear  fractional  LAvy  motions  would 
be  as  the  only  H-ss  si  a-stable  processes  X  of  the  form 

(1.12)  X(t)  =  [  [G(t-u)  -  G(-u)]  dM  (u).  t  €  R. 

J_00  a 

where  G(t  -  •)  -  G(-*)  €  La(IR)  for  each  t  €  1R.  Expressed  in  this  way,  this 
characterization  does  not  involve  the  corresponding  stationary  process  Y,  and 
hence  is  stated  for  all  0  <  a  <  2.  Vervaat  [V2]  derives  such  a 
characterization  under  a  self-similarity  assumption  on  the  kernel  of  (1.12) 
of  the  following  form:  for  all  t,  u  €  R,  c  >  0, 

(1.13)  G(c( t-u))  -  G(-cu)  =  cP{G(t-u)  -  G(-u)} 

and  G(0)  -  0.  Indeed,  (1.13)  implies  Immediately  rhat  G(x)  =  a  +  b  x^. 
This  self-similarity  of  the  kernel  G  also  clearly  implies  the  self-similarity 
of  the  process  X.  Theorem  1  (or  more  precisely,  its  proof  below)  establishes 
this  characterization  when  the  self-similarity  of  C  i«  rep!a<~#»d  bv  that  of  X, 
in  the  case  where  G(x)  =  0  for  x  <  0  or  G(x)  =  0  for 
x  >  0.  We  state  this  characterization  in  the  following. 

Theorem  2.  Fix  0<a<2.  0<H<  1,  H/  1/a.  Let  X  be  a  nondegenerate  H-ss 
si  a-stable  process  of  the  form  (1.12).  If  G  vanishes  on  some  half  line, 
then  X  is  a  one-sided  linear  fractional  L6vy  motion. 
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2.  Proof  of  Theorem  1 . 

Since  the  proofs  are  similar,  we  only  prove  (i)  of  Theorem  1.  Also  for 

notational  simplicity,  we  write  Y  for  Y,  and  M  for  M  . 

La 

00 

Since  Y(t)  =  f  g(t-u)  dM(u)  is  a  nonanticipating  moving  average, 
g(x)  =  0  for  x  <  0.  Since  Jq  |g(x) |a  dx  >  0  (otherwise  Y  =  0  and  thus  X  is 
degenerate),  if  we  define 

tg  =  inf{t  ^  0  |  Leb[x  €  (t.t+6);  g(x)  /  0]  >  0  for  any  6  >  0} , 
then  0  i  t^  <  00 .  We  then  have 

{  CO  ~  00 

Y(t)  =  g(t-u)  dM(u)  =  g(t-(v-t  ))  dM(v-t  ) 

J  —00  J  — on  U  U 


r  <30 

=  |  gt  (t- 
J-€0 


v)  dM(v) . 


where  g  (x)  =  g(x  +  tn)  and  we  have  used  the  fact  that  M  has  stationary 
t0  U 

increments.  By  the  definition  of  t^. 


(2.1) 


Leb(x  €  (0.6);  g  (x)  ^  0}  >  0  for  any  6  >  0, 
0 


and  g  (x)  =0,  x  <  0.  Hereafter  we  write  g  for  g 

0  C0 

From  (1.3)  and  (1.5),  we  have  for  s  <  t. 

,  f  •  rt-u 

(2.2)  X{ t)  -  X(s)  i  [g(t-u)  -  g(s-u)  +  g(v)  dv]  dM(u) 

J -4X)  J  s-u 

=  [  [G(t-u)  -  G(s-u)]  dM(u) . 

J  — €0 


where 


G(x)  =  g(x)  +  [  g(v)  dv,  x  €  R. 
Jn 


Note  that  G(x)  =0,  x  <  0,  since  g(x)  =0,  x  <  0. 

Fix  t  >  0.  Then  we  have  for  any  aj,...,a^,  x^ . x^  €  IR  and  c  >  0, 


-  8  - 


N 

2  an  [X(c(t  -  xn))  -  X(cxn)] 
n=l 


=  f  ]  2  an  CG(c(t  +  xn)_u)  "  G(cxn  ~  U)1 

J-a>  I  n=l 


dM(u) 


■f 


N 

2  an  [G(c(t  +  xn  -  v))  -  G(c(xn  -  v)] 
n=l 


f  dM(cv) 


d  1/a 
=  c 


I. 


N 

2  an  [G(c(t  +  xn  -  v))  -  G(c(xn  -  v))] 
n=l 


dM(v) 


(where  we  have  used  the  1/a-self-simi lari ty  of  M) , 
N 

2  a  [X(t  +  x  )  -  X(x  )] 

,  n  L  v  n7  v  n7J 
n=l 


M. 


N 

2  a  [G(t  +  x  -  v)  -  G(x  -  v)] 
n=l  n  n  n 


dM(v). 


Hence  by  the  self-similarity  of  X  and 

^  co  r  00 

|E[exp{i  f(v)  dM(v)}]|  =  exp{-  |f(v)|a  dv}. 

J  —00  J  —CO 


we  have  for  each  c  >  0, 


I. 


N 

2  an  CG(c( t  +  xn  -  v}}  -  G(c(x  -  v))] 
n=l 


a 


dv 


=  c 


aH 


L 


N 

2  an  [G(t  ♦  x  -  v)  -  G(xn  -  v)] 
n=l 


dv , 


namely,  the  functions  c1/a  H[G(c(t-*))  -  G(c(-*))]  and  G(t-*)  -  G(  —  * )  have 
equal  L^-norms  of  all  linear  combinations  of  their  translates.  Then,  by 
Ranter’s  result  [K],  for  each  fixed  c  >  0, 

(2.3)  c1/a_H[G(c(t-u))  -  G(-cu)]  =  e  [G( t-u-r  )  -  G(-u-r  )]  for  a. a.  u. 

c  c  c 
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where  tc  €  {-1,1},  tc  €  R,  and  they  may  depend  on  c.  We  shall  show  that 

e  =1  and  r  =  0  for  all  c  >  0. 
c  c 

Suppose  rc  >  0.  Then  for  u  satisfying  max{0,t-rc)  <  u  <  t,  (2.3) 


reduces  to 


c1/a_H  G(c( t-u) )  =  0.  a. a.  u, 


since  G(x)  =  0.  x  <  0.  Namely, 


and  thus 


G(x)  =  0  for  a. a.  x  €  {0,  min{ct,  ctc})  =:  Ic> 


g(v)  dv  for  a. a.  x  €  I  . 


The  function  f(x)  :=  J*  g(v)  dv  Is  continuous  and  satisfies  f(x)  =  X*  f(v)  dv 

for  any  x  €  I  .  Thus  f  =  0  on  I  ,  implying  g  =  0  a.e.  on  I  ,  which 
CO  c 

contradicts  (2.1). 

When  we  suppose  <  0,  it  is  enough  to  take  u  satisfying  max{t,-Tc)  <  u 
<  t  -  in  order  to  get  a  contradiction.  Namely,  for  such  u,  we  have 
G(t-u-Tc)  =  0,  so  that  G(x)  =  0  for  a. a.  x  €  (0,  min{t,  — r  }).  Thus  =  0 


for  any  c  >  0,  and  we  obtain 

cl/a  ^[G(c(t-u))  _  -  e^[G{t-u)  -  G(~u)]  for  a. a.  u. 

We  thus  have  for  a. a.  u  €  (0,t),  c^°  ^  G(c(t~u))  =  G(t-u),  namely 

(2.4)  G(cx)  =  &c  cH  G(x)  for  a. a.  x  €  (0,t). 

Choose  0  <  a  <  b  <  t  and  integrate  (2.4);  then 


f1*  .  H-l/a  f* 

G(cx)  dx  =  e  c 


G(x)  dx 


so  that 


(2.5) 


f  G(x)  dx  =  ec  cH_1/a+1  I**  G(x)  dx. 

J  t-s  r*  J  n 


thus  have  the  same  contradiction  as  before  for  each  c  >  0.  If  G(x)  dx  *  0 

& 

for  some  0  <  a  <  b  <  t.  then  we  see  from  (2.5)  that  e.  is  a  continuous 

c 

function  of  c  >  0  and  thus  t  =  t,  =  +1  for  all  c  >0. 

c  1 

Therefore,  we  obtain 

(2  6)  G(cx)  =  cH  G(x)  for  all  c  >  0  and  a. a.  x  €  (O.t). 

By  the  same  reasoning  as  before,  we  see  that  G(Xq)  ^  0  for  some  Xq  >  0 
independent  of  c  >  0,  We  thus  have  from  (2.6), 

(2.7)  G(x)  =  (x/Xq)**  1/a  G(Xq)  =  a  x^  1/a,  for  all  x  >  0. 

where  a  =  G(Xq)  ^  ^  0.  We  know  that  G(x)  =  0  for  all  x  <  0.  and 

assuming  G(0)  =  0  does  not  change  the  distribution  of  X.  So,  we  have 

(2.8)  G(x)  =  a  x^  for  all  x  €  IR. 

Noting  that  any  H-ss  si  process  with  H  >  0  satisfies  X(0)  =  0  a.s.  (see, 
e.g..  [VI]).  we  conclude  from  (2.2)  and  (2.8)  that 

X(t)  =  a  [  [(t-u)“'1/a  -  (-u)“_1/a]  dM(u) , 

+  -40 


which  completes  the  proof. 
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